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should  read:  "  R  ^  ^  ^  S  ^ 

C  U  u  dt  u 

"  atmosphereic  "  should  read;  "  atmospheric  " 
"  equation  (12)  "  should  read:  "  equation  (18)  " 

"  equation  (1)  "  should  read:  "  equation  (19)  " 
should  read:  "  s  =  Ae^^*  +  Be  " 


should  read:  "  where  A  and  B  are  arbitrary  constants  and:  " 
should  read:  "  X  =  Ik  +  Ik 


I  =  w  +  i .« 


"  e-’"*  or  e*'*^*  " 


m 

should  read: 


iXt 


or  e 


-iXt 


"  in  Equation  (35)  "  should  read:  "  in  Equation  (31)  " 
Equation  (59)  should  read:  "  s^  +  ^ 


f  j(fi ) 


f^{B) 


s  0  " 
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"  Figure  10  "  should  read:  "  Figure  30  " 
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"  Equation  (82)  "  should  read:  "  Equation  (84)  " 
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AB3TRACT  ^ 

XWptpan  pHaaiiH  lolutiont  of  Reynold!  equjLtlon  for  cylindrical 
Journal  bearlngi  in  which  the  velocity  of  the  Journal  center  le  considered; 
however,  the  inertia  terms  are  neglected.  Solutions  are  obtained  for  balh 
compressible  and  incompressible  fluids.  Analysis  of  the  steady  state 
forces  that  exist  in  compressible  fluid,  externally*pressuriEed  cylindrical 
Journal  bearings  with  one  plane  of  feeding  is  also  included.  From  these 
analyses  the  gradients  of  the  radial  and  tangential  forces  with  respect  to 
displacement  and  velocity  are  obtained.  The  equations  of  motion  for  small 
and  large  oscillation  In  terms  of  these  forces  and  gradients  are  set  up.  and 
stability  criteria  are  established  for  small  oscillations  and  also  for  large 
oscillations  of  a  noassless  rotor.  Sample  calculations  for  threshold  of 
instability  are  given  and,  where  existing  experimental  data  is  available, 
theoretical  predictions  are  compared  with  results  attained  in  practice. 
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.WYRODUCTION 


Tha  dynamic  charaetarlatlea  of  rotating  aquipmant  ara  baeomiag  vary 
much  mora  Important  with  tha  praaant  trend  towards  higher  apaads  and 
praeiaion  equipment  auch  aa  turbinaa.  compraaaors,  gyroai  grinders,  ate. 
Numerous  examples  can  be  cited  where  equipment  was  unacceptable  or 
actually  failed  due  to  instability,  excessive  vibrations,  operating  speed  co¬ 
inciding  with  system  criticals,  and  other  dynamic  malfunctions. 

In  recent  years  it  has  been  pointed  out  that  bearings  have  major  influence 
on  dynamic  raaponse  of  rotors.  They  possess  both  atiffneaa  and  damping 
and  thus  govern  the  system  criticals  and  amplitudes  of  vibration.  It  has 
also  bean  shown  that  forces  generated  in  the  fluid  films  of  the  bearings  may 
excite  lateral  vibration  in  a  rotor  and  cause  it  to  become  imstabla.  This 
instability  may  be  so  violent  that  the  bearings  fail  within  minutes  and  often 
damage  the  whole  machine.  The  importance  of  this  problem,  therefore,  has 
attracted  many  investigators  since  it  was  first  reported  by  Newkirk  and 
Taylor  in  1925.  However,  tha  subject  is  very  complicated  and  it  is  not  sur¬ 
prising  that  many  observations  and  even  the  terminology  ara  controversial. 

A  brief  history  may  be  in  order  so  as  to  highlight  soma  of  these 
observations.  Newkirk  and  Taylor  (1)  reported  a  violent  whip  of  shafts 
which  was  caused  by  the  oil  film  in  the  bearing.  They  found  that  the  whip 
started  at  speaditwlce  the  critical  and  persisted  at  higher  speeds.  Often 
the  instability  could  be  initiated  merely  by  shock.  Robertson  (2)  made 
certain  assumptions  and  reached  the  conclusion  that  the  rotor  is  inherently 
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uniUble  at  all  apaada  of  rotation.  Poritaky  (3)  ahowad,  by  Introducing  an 
additional  force  to  the  oil  film  forcoi  that  the  Inatablllty  oecura  only  at  apaada 
axcaeding  t^vlca  the  flrat  critical  apaad.  Soaker  and  Sternlicht  (4)  pointed 
out  that  the  additional  force  introduced  by  Poritaky  (3)  in  reality  doea  exiat 
and  it  correaponda  to  the  t:<dial  component  of  the  hydrodynamic  force.  Ref.  4 
further  pointa  out  that  the  threahold  of  Inatablllty  la  dependent  on  the  magnitude 
of  this  force  and  that  under  aome  conditiona  inatablllty  may  atart  even  at 
£ero  apaad.  It  alao  aho'va  that  the  analyala  appliaa  equally  to  compraaaible 
and  incompreaaible  hydruu;  .lamic  liibricatlon.  provided  the  hydrodynamic 
force  a  are  properly  calculated  for  the  two  caaea.  Thla  important  obaer\'atlon 
waa  verified  experimentally  in  Raf.  (4)  and  more  recently  Fiaeher,  Cherubim 
and  Decker  (5)  have  confirmed  thla  in  theii  experimenta.  On  the  other  hand. 
Newkirk  and  Lewie  (6,  7)  and  Jukla  Hori  (8)  reported  experimental  caaea  in 
which  the  rotating  apeed  reached  five  or  aix  timea  the  critical  before  the 
Inatabllity  developed. 

Newkirk  and  Taylor  (1 ,  7)  reported  alao  that  the  violence  c(  whip 
increaaea  with  apeeda  higher  than  twice  the  critical  ■  Hori  (8)  and  Siemllcht 
(9)  obaerved  that  in  aome  caaea  the  amplitude  of  whip  decreased  ae  the  speed 
Increased.  Pinkua  (10)  found  even  more  complicated  examplea  in  which 
whip  disappeared  and  appeared  again  as  speed  was  raised. 

The  effect  of  fluid  \riacosity.  turbine  or  motor  electrical  loads  on 
instability  are  even  more  controversial. 


Another  claei  of  inttabillties.  often  referred  to  &e  "Fluid  Criticale"i 
are  rotor  initabilitlee  of  limited  speed  range  due  to  fluid-film  action.  Theie 
have  first  been  discovered  by  Stodola  (11)  sind  further  studied  by  his  pupil 
Hummel  (12).  Cameron  (13)  pointed  out  that  aside  from  the  two  critlcals 
discussed  by  Hummtli  there  also  Is  a  third  critical  in  the  region  of 
eccentricity  ratio  less  than  0.7. 

These  observations  point  out  some  of  the  complexities  and  controversies 
in  the  subject.  In  this  study  the  authors  give  a  more  comprehensive  theory 
of  stability  which  explains  some  of  the  inconsistencies  and  gives  a  better 
insight  into  the  subject  matter. 

(In  order  to  clarify  the  terminology  employed  in  this  paper,  a  section 
entitled  "The  Mechanics  of  Hydrodynamic  Instability"  is  given  which  defines 
some  of  the  terms. ) 

The  report  gives  a  number  of  solutions  of  Reynolds  equation  for 
isothermal  compressible  and  Incompressible  lubrication  of  plain  cylindrical 
Journal  bearings.  The  analysis  neglects  inertia  terms  but  considers 
velocity  of  the  journal  center.  In  the  case  of  compressible  lubrication, 
p  3  14.7  psia  is  assumed  at  e  e  t  L/2.  Few  cases  are  also  given  for  120° 
partial  Journal  bearing  operating  with  air  as  lubricant.  Results  are  given 
in  dimensionless  foim\  for  the  radial  and  tangential  forces  as  functions  of 
(E.  L/D). 

Several  cases  are  also  solved  for  the  compressible  externally 
pressurized  cylindrical  Journal  bearing  with  central  single  plane  of  feeding, 
operating  under  non-rotating  steady  state  condition.  Dimensionless  force 
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is  given  as  a  function  of  (t,  L/D,  p  ,  p  , 

sat  opt 

A  brief  discussion  is  next  included  to  show  that  instability  is  induced  by 
forces  nomaa.1  to. the  radlll  displacement  of- the  Journal.  Following  this,  equations 
and  methods  of  calculation  are  given  for  the  threshold  of  half  frequency  whirl 
for  both  vertical  and  horisontal  symmetrical  rotors-  Equations  are  also  given 
for  calculating  the  amplitudes  of  forced  vibration  of  vertical  shafts.  It  is 
fur  ler  shown  ho\  le  fluid  film  and  rotor  resonance  in  rigid  vertical  shafts 
can  be  calculated. 

Equations  of  motion  are  set  up  for  small  oscillations  about  a  position  of 
equilibrium.  The  equations  are  in  terms  of  displacements  and  the  hydrodynamic 
radial  and  tangential  forces.  Solutions  of  these  equations  exist  in  which  time 
enters  as  an  exponential.  Several  solutions  are  given  for  a  symmetrical  rotor 
supported  by  ^  elastic  shaft  and  fluid  film.  Results  give  the  threshold  of 
resonant  whip  as  a  function  of  eccentricity  ratio  and  dimensionless  parameter  E 

The  more  general  equations  for  large  displacements  are  also  set  up  in 
terms  of  eccentricity  ratio  and  attitude  angle  These  equations  form  a  set  of 
two  non-linear  differential  equations,  each  of  order  3.  In  the  case  of  small 
rotor  mass,  the  terns  involving  m  can  be  neglected  and  the  differential  equations 
reduce  to  a  second  order  non-linear  system.  The  integration  of  these  can  be 
carried  out  by  a  graphical  method  and  the  steps  are  outlined. 

For  several  cases  where  ('xperimental  data  was  available,  comparison 
between  theory  and  practice  is  given.  Explanation  is  given  for  several  of  the 
reported  observations  which  appear  to  be  irreconcilable  It  is  shown  that  in 
reality  they  provide  consistent  verification  of  the  theory  and  thus  enlarge  our 
available  knowledge  for  a  bener  understanding  of  Journal  bearing  instability. 
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U  THE  MECHANICS  OT  HYDRODYNAMIC  INSTABILITY 

KydrodyTLvnle  IniUbllity  It  c&uted  by  foreet  gtntrattd  In  tb«  fluid  film 
of  the  bearing)  to  directed  with  reepaet  to  the  ehaft  dleplaeement  at  to  propel 
it  In  Iti  whirling  motion.  To  vieualiae  better  the  forces  which  cause  this 
whirling  motion,  consider  Fig.  1  which  shows  the  equilibrium  attitude  of  a 
shaft  rotating  in  a  bearing  and  supporting  a  load  W.  As  this  It  the  position  of 
equilibrium,  the  resultant  force  of  the  fluid  film  on  the  shaft  Is  the  force  F, 
which  must  be  equal  and  opposite  to  W.  The  important  thing  to  notice  in 
Fig.  1  is  that  the  force  F  Is  not  in  the  direction  of  centers  O'O  but  at  an 
angle  ^  to  it. 

Now  let  us  consider  a  rotating  shaft  which  carries  no  load,  (W  b  o,) 
and  suppose  that  It  is  momentarily  displaced  from  0  to  O',  as  shown  in 
Fig.  2.  The  film  under  these  conditions  sicerts  a  resultant  force  F  on  the 
shaft,  Just  as  it  did  in  Fig.  1.  However,  in  this  case  tiiere  Is  no  opposing 
force,  W,  so  the  force  F  must  be  spent  in  accelerating  the  shaft  and  over¬ 
coming  the  frictional  drag  of  its  resulting  motion.  The  movement  of  the 
shaft  center  O'  in  response  to  the  force  F  obviously  will  not  be  toward  the 
shaft  center  O-  Rather,  the  shaft  center  will  be  forced  to  progress  in  an 
orbit  around  the  bearing  center  and  as  long  as  the  centers  do  not  coincide, 
some  force  F  will  be  generated  by  the  rotation  of  the  shaft  and  the  whirl  will 
continue.  Whether  the  whirling  motion  becomes  more  pronounced,  continues 
at  the  same  amplitude,  or  dies  out.  depends  upon  the  angle  4*  and  the  damping 
characteristics  of  the  bearing  and  shaft  system. 
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The  whirl  frequency  ii  eet  by  the  speed  et  which  the  theft  etn  pump  the 
fluid  around  in  the  clearance  and  maintain  the  preieure  pattern  which  produces 
the  driving  force  F.  Aeeuming  laminar  conditions  and  neglecting  flow  due  to 
any  pressure  gradients  (which  In  a  lightly  loaded  bearing  are  very  small),  the 
average  lubricant  velocity  ir  half  the  peripheral  speed  and  in  the  direction  of 
shaft  rotation.  Hence,  the  greatest  frequency  at  which  a  pressure  pattern 
can  progress  around  the  bearing  is  half  the  shaft  speed.  This  is  the  speed 
and  direction  at  which  the  shaft  vibration  occurs.  We  shall  refer  to  this 
kind  of  instability  as  "Half-Frequency  Whirl".. 

Finally,  consider  the  case  illustrated  in  Fig.  ?.  in  wb*ch  the  shaft 
supports  a  load  W  in  the  equilibrium  position  O'.  Suppose  the  shaft,  by  some 
external  shock,  is  momentarily  given  a  secondary  displacement  to  a  new 
position  O".  The  fluid  force  F  corresponding  to  this  new  position  of  the 
shaft  is  no  longer  equal  and  opposite  to  W.  The  vector  difference  between 
the  forces  F  and  W  is  a  force  F^  which  can  cause  the  shaft  center  to  whirl 
(in  the  direction  of  shaft  rotation)  around  its  equilibrium  position  O'  at  a 
speed  nearly  equal  to  half  the  shaft  rotational  speed  •  Just  as  in  the  case  of 
the  unloaded  shaft  discussed  above.  As  before,  the  persistence  of  such  a 
whirling  motion  will  depet  d  upon  the  damping  characteristics  of  the  system. 

Another  phenomenor,  system  resonance,  can  Join  with  half  frequency 
whirl  to  produce  a  vigorous  vibration.  When  the  rotational  speed  is  about 
twice  the  actual  system  fir>it  critical,  the  system  will  build  up  in  resonance 
at  a  frequency  equal  to  the  system  first  critical  frequency.  This  form  of 
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reionanca  li  raferrad  to  aa  "Raaonant  Whip"  and  may  ba  dafinad  aa  a  raaonant 
vibration  of  a  Journal  In  a  fluid  film  baaring  whicht  for  lovp  aceantrielty  ratloi, 
lata  in  at  approximately  twice  the  actual  firat  ayatam  critical  and  peralata 
at  higher  apaada  with  the  frequency  of  vibration  approximately  equal  to  the 
firat  ayatam  critical  regardlaaa  of  running  apead.  Hera  too  the  motion  of 
the  ahaft  canter  la  in  the  aama  direction  aa  ahaft  rotation.  Raaonant  whip 
ia  a  aalf-aupportad  vlbrationi  aa  ia  half  frequency  whirl.  In  the  caaa  of 
raaonant  whipi  the  vibration  la  aupportad  by  the  fluid  film  action  while  the 
frequency  la  controlled  by  the  ayatam  critical  apead. 
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UI  CALCULATION  OF  FLUID  TILJ^  FORCES  IN  SELF  ACTING  BEARINGS 
The  general  form  of  the  Reynolds  equation  may  be  written 

It  ^  It  p||  ‘  ^  ^  +  2  p  V  +  2h  |£  (1) 

In  thir/  equi;'  .on  U  is  the  relative  tangential  velocity  between  the  Jouvnal  and 
heaving  surfaces  and  V  is  the  relative  normal  velocity  between  them. 

Equation  (1)  then  contains  the  terms  which  determine  the  instantaneous 
velocity  of  toe  Journal  center.  The  fluid  inertia  terms i  however,  are 
neglected  in  the  derivation  of  Equation  (1),  In  accordance  with  past  analytical 
and  experimental  evidence  which  show  their  effect  to  be  small,  at  least  at 
mod/  rate  speeds. 

The  hydrodynamic  pressures  (and  hence,  by  integration,  the  fluid 
film  forces)  are  obtained  from  the  Integration  of  Equation  (1)  using  the 
appropriate  boundary  conditions. 

For  a  full  circular  bearing,  Equation  (1)  is  subject  to  the  boundary 
condition  of  ambient  pressures  at  the  edges.  Thus: 

p  ( 0,  L/2)  =  p  (0,  -L/2)  ■  p  (2) 

a 

We  will  now  further  examine  Eqttatlon  (1),  first  for  the  case  of 
an  incompressible  fluid,  then  for  the  case  of  a  compressible  fluid  under 


isothermal  conditions. 
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A.  Incotnpretsible  Fluid 

If  the  fluid  !■  ineompretelble»  the  density  is  assumed  constsnt  and  can 
therefore  be  faotored  out. of  Equation  (1),  which  then  reduces  to; 


d 


(1  +  E  cosei^ 


+  R' 


(1  +  e  cosS)* 


‘Ml 


e(w-2a) sin  6  *  Zt  cos  & 


Also,  in  the  case  of  an  incompressible  fluid,  cavitation  may  occur.  If  the 
pressure  "p"  *ends  to  take  on  values  below  ambient  over  an  area  A 
extending  to  the  ends  of  the  bearing,  then  cavitation  occurs  over  A^.  Equation 
(3)  then  applies  outside  A^.  along  with  the  additional  boundary  conditions: 

pop  over  A 


(4) 

^  »  0  over  Ci  the  boundary  of  A^ 

(where  n  Is  the  directional  co-ordinate  norrna!  to  C) 

If  cavitation  were  to  occur  over  an  Internal  area  not  extending  to  the 
bearing  ends,  pressure  below  ambient  could  occur  and  the  onset  of  cavitation 
would  be  governed  by  the  boundary  conditions: 
p  =  0  over  A 


0  over  C 


{4a) 


It  turns  out,  however,  that  no  such  situation  occurs  and  Equation  (4) 
applies. 

For  incumpresslble  fluids,  p  in  Equation  (1)  may  be  used  to  denote  the 
pressure  above  aunblent  (gage  pressure)  in  which  case  we  set  p  (gage)  *  0 

A 
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Thfi  integration  of  Equation  (3)  yields  the  pressure  distribution  in  the 
bearing  corresfonding  to  the  conditions  at  which  the  integration  is  performed. 
By  a  further  integration  of  the  pressures  outside  of  Q,  one  obtains  the  net 
resultam  force  which  the  fluid  film  exerts  on  the  Journail  and  thus,  one  may 
evaluate  the  components  and  F^.  of  this  force  in  the  radial  and  tangential  : 
directions  respectively.  ' 

From  dimensional  analysis  considerations,  these  •components  can  be 
expressed  in  the  following  form: 

V  V. 


F^  =  -uLU 


S  f  (c  ~  —I 

‘r  U  ■  U'  d' 


(S) 


ulLU  §  *  r  (e, 


r  t  L. 

u  ’  u’  d' 


where  {-  Ce)  and  (a  Cea)  are  the  radial  and  tangential  components, 
respectively,  of  the  Journal  center  velocity,  and  where  and  are 
dimensionless  and  depend  on  the  dimensionless  variables  shown. 

Note  that  the  dimensionless  quantities  V^/U  and  V^/U  may  be  replaced 


as  follows: 


rIi 

C  '  U 


.  ide 

u  dt 


£ 

w 


C‘  U 


1  ^ 

u  dt 


a 

u 


For  any  given  L/D  ratio  then,  Equation  (5)  may  be  written: 


F 

r 


■uLU  ^  •  £  (e.  e'.  o') 


^LU 


#  *  f  it.  o') 


{5a) 
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The  negative  sign  in  the  equation  ia  due  to  the  fact  that  ia 


measured  poeitive  in  the  direction  of  increasing  radial  displacement,  while 
(at  least  for  c  0)  the  radial  force  turns  out  to  point  in  the  direction  opposing 
the  radial  displacement. 

The  right  hand  member  of  Equation  (3)  is  linear  in  w.  a.  i.  It  would 
appear  from  this  that  ReyTiolda  equation  for  incompressible  fluids  can  be 
integrated  separately  first  for  a  statioMry  Journal  center,  (u  /  0,  I  b  a  »  0), 
then  for  a  radial  velocity  of  the  Journal  center  (u)  >  d  «  0,  E  U),  and  finally 
for  a  tangential  velocity  of  the  Journal  center  (u  s  e  >  0.  d  /  0).  and  hence 
by  superposition,  the  pressures  (and  thus  the  forces  F^  and  F^}  obtained  for 
the  joint  effects  of  U,  and  V^.  This  condition,  however,  is  vitiated  by 
the  non-linear  requirement  of  the  boundary  condition  given  in  Equation  (4) 
and  the  fact  that  the  curve  C  will,  in  general,  change  in  a  manner  depending 
on  all  three  velocities  U,  and  V^. 

Nevertheless,  examination  of  the  right  hand  side  of  Equation  (3)  shows 

that  at  least  U  and  enter  in  the  form  (eCU  -  2RV^)  sin  6,  or  its  equivalent 

form,  (w  -  2d)  sind.  Suppose  that  Equation  (l)iB  integrated  for  a  given 

L/0  and  the  forces  F  and  F^  evaluated  for  >  0  yielding 

r  t  t 


-X  w  f^  (e,  £') 
X  u  f ^  (e,  e’) 


(6) 


where  X  u  has  the  dimensions  of  a  force,  x  being  defined  by: 


X  » 


uLR 


(7) 
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Then,  for  non-zero  or  a.  the  forces  may  be  derived,  without  further 
integration,  by  replacing  w  by  w  (1  -  2a*).  At  the  same  time,  the  dimensionless 
parameter  e'  must  be  divided  by  (I  -  2a').  Hence,  for  i  0,  the  fluid  film 
forces  become: 


»  -X  «  (1  -  2a»)  f^  {ei  e') 
and  »  X  w  (1  -  2q')  f^  (€.  f) 

whe  re 


(8) 


(8a) 


For  the  case  of  a  full  circular  bearing  and  incompressible  fluids, 
discussed  above,  numerical  integrations  of  the  Reynolds  equation  have  been 
carried  out  for  L/D  ratios  of  0.  5,  i .  0  and  1 . 5,  and  for  a  wide  range  of 
Journal  center  speeds.  The  resulting  fluid  film  forces  are  given  In  Tables  1 
and  la,  the  symbols  in  these  tables  corresponding  to  the  nomenclature  given 
in  Figure  4. 

For  the  case  of  1./D  •  1,  these  forces  and  attitude  angles  have  also 
been  plotted  in  Figures  5  and  6,  on  the  eccentricity  ratio-journal  center 
velocity  plane  Note  that  these  curves  represent  an  extension  of  the  conventional 
one  parameter  steady-state  solution  to  the  two  parameter  dynamic  condition. 

The  conventional  steady-state  solution  appears  in  this  more  general  repre¬ 
sentation  as  the  particular  case  where  e'  «  0. 


i- 
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B.  Compreiiibl«  Fluldi  -  Iiotherraal  Condltiofti 

When  the  bearing  fluid  it  a  gae  and  the  ratio  of  the  film  to  ambient 
pref^iuree  le  largei  compreeiibility  effects  may  no  longer  be  neglected  and 
the  more  general  form  of  the  Reynolds  equation  given  In  Equation  (1)  applies. 
The  right  hand  side  of  Equation  (1)  includes  the  non-steady  flow  term: 


This  term  may  be  of  importance  in  cases  of  high  frequency  oscillations 
of  the  Journal  center.  At  low  frequency  oscillations  of  the  Journal  center, 
however,  it  is  of  smaller  magnitude  than  the  other  terms  on  the  right  hand 
side  of  the  equation.  Furthermore,  its  inclusion  in  the  analysis  would  greatly 
complicate  the  work  by  introducing  a  third  parameter  (time)  In  the  solutions. 

In  the  present  analysis  for  compressible  fluids,  the  non-steady  flow  term 
has  been  neglected.  More  rigorous  examination  of  its  effect  on  bearing 
stability  Is  planned  for  the  future. 

We  have  also  assumed  in  the  present  analysis  that  Isothermal  conditions 
prevail  in  the  flxiid  film  and  this  conforms  with  the  results  of  numerous 
experiments  (set.  for  example.  References  S.  15  and  lb). 

Under  isothermal  conditions,  the  equation  of  state  is: 


p  ■  ft  T  p 


(9) 


where  ft  is  the  gas  constant  and  T,  the  absolute  temperature. 

Eliminating  p  from  Equation  (1)  by  means  of  Equation  (9)  and 
dropping  out  the  non-steady  flow  term,  the  Reynolds  equation  becomes: 


+  2  p  V 


(10) 


>14. 


This  equation  it  non-linear.  Thus,  the  method  of  tuperpooition  (which 
was  proved  to  be  Invalid  'in  the  incompressible  cate  on  account  of  the  non- 
linear  boundary  conditiont  Equation  (4)^  is  Invalid  here  too  because  of  the 
non-linearity  of  Equation  (10). 

The  boundary  condition  of  atmosphereic  pressure  at  the  bearing  ends, 
Equation  (2),  still  holds,  however,  absolute  pressure  in  the  bearing  region 
may  now  became  sub  .atmospheric  (though,  of  course,  never  negative). 

Once  again,  as  in  the  incompressible  case,  the  fluid  film  forces  may  be 
expressed  in  terms  of  their  radial  and  tangential  components  and  F^, 
respectively.  Here,  of  course,  the  dimensionless  functions  f^  and  f^  have 
to  be  calculated  anew,  this  time  by  integration  of  the  non-linear  Reynolds 
equation  (Equation  10). 

Numerical  integrations  have  been  carried  out  for  full  circular  bearings 
(with  compressible,  isothermal  flow)  for  L/D  ratios  of  0. 5,  1.0  and  1 .  5. 
and  a  wide  range  of  Journal  center  velocities.  These  are  given  in  Table  11. 
For  purposes  of  comparison,  similar  additional  cases  have  been  calculated 
for  a  120  deg.  partial  bearing  with  an  L/D  ratio  of  0.  3  and  these  are  given 


in  Table  in. 
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IV  CALCUtATION  OF  rORCKa  ON  JOURNAL  IN  EXTEK^^ALLY- 

PRESSURIZED  BKARINOS 

Anjilyi.ei  have  been  made  based  on  (Ref.  15,  16)  for  the  resultant 
hydrostatic  force  acting  on  a  Journal  in  a  cylindrical  bearing  with  central 
plane  of  feeding.  Results  are  given  in  terms  of  L/D  ratio,  supply  pressure 
p^,  ambient  pressure  p^  Sind  Z^. 

The  governing  equation  which  applies  for  an  incompressible  fluid 
non- rotating  journal  and  small  eccentricity  ratio  is: 


^  3  Sh 

*  cs 


a  0 


Since  ^  -^0  as  e — $0,  for  small  e.  of  the  same  order  of  magnitude 
te  e<  This  means  that  ^  ^  ^  is  of  the  order  of  (e)*.  Thus  for  e  <<  1, 
equation  (11)  reduces  to 


a  0 


Ec<  1 


(lla) 


In  the  case  of  eompreseible  fluids  under  Isothermal  flow  conditions 


equation  (11)  la  of  the  form 


ph* 


1  9  ' 

+  h*  ^  ! 


57  [P  57, 


B  0 


Performing  the  indicated  operations,  the  equation  expands  to 


The  film  thickness  equation  may  be  written  at 


a  0 


(12a) 


h  B  C  (1  -  E  cos  6] 
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Thli  equation  ie  of  the  iiune  fom  at  Equation  (iZa) 

Applying  the  aeparatlon  technique  to  Equation  (12d)  yiclde  two  ordinary 
differential  equation!  which  can  be  eolyed  for  ^  a). 

It  ehould  be  pointed  out  that  for  rigoroui  analueie  of  the  extemally- 
preeeuriced  bearing,  the  equation  of  continuity  of  flow  between  the  feeders 
and  the  bearing  must  also  be  aatiefied.  Since  the  feeders  and  bearing 
clearance  are  restrictors  in  series,  they  Jointly  affect  the  bearing  per¬ 
formance.  The  calculations  presented  here  take  account  of  this  and, 
therefore,  the  results  are  in  terms  of  the  dimensionless  fe-iding  parameter 
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In  general  one  le  intereeted  in  getting  maximum  etiffoeea  from  a  hydro* 

•tatlc  bearing.  Thle  le  aecompllahed  by  optlmialng  the  parameter  veraui 

the  dlmenilonleee  eprliig  c<>netant  1.  Oh  the  left  hand  tide  of  the  eptlmom 

tee  much  reetrlction  li  preeent  in  the  feeding  eyetemi  while  tiie  reverie  le 

true  on  the  right  hand  elde  of  Z^ 

t  opt 

The  procedure  for  ealeulatiene  la  outlined  in  order  to  clarify  the  eubjeet 
matter; 

1 )  Olven  L/D,  p^  and 
2}  Aeaume  valuee  of  Z^ 

3)  Calculate 


4)  rind  g  (Y,  L/D) 


where  i/  and  ^  are  the  error  functlone 
*  (*)  “  /*  •* 

f  /*  .-‘’d, 

1  O 

5)  Calculate  dlmenalonleae  etiffneae  k 


(15) 
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*nd  optimise  k  by  vsrylnj  2^  (•««  Figure  7\. 
Stifiixesi  can  be  calcuUt^d  from 


£p  ^ 
C 


(16) 


6)  C&lculste  dixnena ionises  pressure,  p  vs  e 


Q 


(17) 


The  results  sre  given  In  Fig.  8-iO  ior  three  dl^e rent  supply  pressures 
and  three  L/D  rntlos. 


7)  Restrictor  ealeulntions 


(18) 


where 

■  dlsehArfe  eoelfleUnt  (iir.0.6) 

d  ■  orifice  diemeter 
n  B  number  of  fesders/beerlnf 
m  «  number  of  orifice  s/feeder 


Several  important  conclusions  can  be  drawn  from  this  analysis. 

a)  For  a  non-routing  estemaUy  preesuHasd  Journal  bearing,  the  dlspUeement 
is  in  the  direction  of  load  (attitude  angle  is  sero). 

b)  The  displacement  Is  proportional  to  lead  for  low  eccentricity  ratios 

0.  5.  Thus,  the  fluid  stiffness  is  a  constant.  This  leads  to  the 
conclusion  that  for  a  given  mass  rotor  the  critical  speed  of  the  system 
Is  independent  of  eccentricity  ratios  for  e  c  0.5. 
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c)  The  critical  speed,  can  be  raised  by  Increailng  the  supply  preeiure. 

This  is  confirmed  by  our  experiments  and  those  reported  in-Ref.  5. 

d)  From  the  standpoint  of  maximum  load«  it  appears  that  there  is  an 
optimum  L/D  ratio.  This  shows  up  clearly  in  Figure  9*  where  It  is 
seen  that  (for  the  values  of  the  parameters  used  in  that  figure) 

a  1.  Thusi  Figure  9  shows  that  a  bearing  with  L/Db  z 
which  has  twice  the  area  of  the  optimum  bearingi  carries  only 
of  its  load.  Similarly!  a  bearing  with  .L/D  a  i/z,  which  has  half 
the  area  of  the  optimum  bearing,  carries  of  its  load. 

e)  Once  the  optimum  is  obtained  one  has  still  ire»  >/m  of  choice  of  . 
several  parameters  in  equation  (12).  e.g..  d.  C.  n  a. ' m. 

It  has  been  shown  by  several  invsstlgavrs  that  the  hydrodynamic  effects 
are  very  small  when  siqierimpossd  on  the  bydirostatle  effects.  This  is 
especially  true  when  the  hydrostatic  supply  pressure  is  high  (p^/p^>  3). 

There  is.  however,  a  small  attitude  angle  when  rotation  is  coxisldered.  The 
attitude  angle  decreeees  with  Increaee  in  eupply  pressure  for  under  this 
condition  hydrostatic  effects  predominate.  Experimental  reeulte  confirming 
this  point  are  shown  in  Fig.  1 1  ■ 

From  Fig  8-11  one  can  calculate  the  optimum  force  and  stiffness  as  a 
function  of  e,  L/D.  p  /p  -  From  Fig.  1  l  one  can  estimate  the  attitude  angle  as 

9  A 

a  function  of  p^  and  in  turn  obtain  the  radial  and  tangential  force  and  stiffness 
for  the  hydrostatic  bearing  in  which  hydrodynamic  effects  due  to  rotation  are 
considered- 
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V  INSTABOJTY  INDUCED  BY  FORCES  NQRX4AL  TO  RADIAL  DISPLACEMENT 
In  the  preceding  eeetlone  we  have  obtained  the  fluid  film  forcee  that  act 
on  the  Journal.  Before  proceeding  to  uie  theae  In  the  atabillty  inveitigationa. 
it  la  of  Intereat  to  point  out  the  effect  of  the  force  component  normal  to  the 
radial  displacement. 

The  equAtlona  of  motion  of  a  balanced  rotor  with  a  restraining  force  and 
a  force  normal  to  the  displacement  from  equilibrium  and  proportional  to  It 
are: 

mx  'll  -kx  +  k.y 

;  O’) 

my  *  -ky  -  k^x 

Let: 


z  «  X  +  iy  (where  1  *  5  (20) 

Multiplying  the  second  equation  (1)  and  adding  It  to  the  first  onet  there 
results: 


me  *  •  (k  +  Ik^  )  e 

whose  solutions  are: 


e  «  Ae  +Be 

where  A  and  B  are  arbitrary  constants.  J  Indicates  complex  components,  and 


^  ^2 
m 


w  +  1  w. 


Since  >  is  a  complex  number,  it  follows  that  either  e*^^^  or  will  become 
infinite  with  time.  Hence,  a  rotor  with  a  force  normal  to  Its  displacement 
is  Inherently  unstable. 
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VI  THBXSHOt^  FOR  HALT  rR£QPENCY  WfORL 

Employing  the  hydrodynemic  forcei  generated  within  the  fluid  filzni  it 
haa  been  ahown  (Ref.  3  and  4)  that  the  thre ahold  of  inatability  may  be 
repreaented  by 

mM*{^  +  ^)<4  (21) 

2 

where  I/k^  i*  the  radial  fluid  film  reailianee.  la  Ref.  4,  it  waa  ahown  that 
thia  inequality  appliea  to  bearinga  which  poaaeaa  a  high  degree  of  aymmetry. 

It  waa  alao  ahown  in  that  reference  that  the  reaulta  apply  to  both  compreaalble 
and  Incomprecaible  fluida.  Thia  point  waa  further  confirmed  es^rimentally 
in  Ref.  5. 

In  Ref.  2  it  la  ahown  that  neglecting  friction  foreaa.  the  reaultant  film 
force  due  to  Journal  rotation  in  a  plain  bearing  la  at  right  anglea  to  the 
eceentriclty*  and  ia  Juat  balanced  by  reaultaid  force  due  to  whirling  when  the 
whirling  frequency  ia  ene«half  rotational  frequency.  Thia  reault  repreaanta 
the  i^iper  limit  of  whlrlUig  frequency  for  an  unloaded  ideal  bearing.  The 
■ame  reault  can  be  deduced  more  directly  by  conaldering  the  continuity 
condition  that  the  volume  of  lubricant  paaaing  through  the  film  at  aoma  point 
A,  Fig.  12  muat  equal  the  volume  paaaing  point  B  plua  the  volume  required 
to  fill  the  void  of  the  receding  Journal.  In  the  limiting  condition  no  preaeure 
can  exlat  in  the  filmi  he*'xe  aide  leakage  ia  aero,  and  it  followa  that 

^(C+e)  1  ^(C-e)  +  2ROe 

Hence: 

O  a  u/2  (22) 

The  aame  reault  ia  obtained  regardleaa  of  where  pointe  A  and  B  are  taken. 
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vn  FORCED  VIBRATION  OF  VERTICAL  ROTOR  {SIMPLin3SD  ANALYSIS) 
Let  ue  consider  a  vertical  rotor  supported  by  two  plain  Journal  bearings 
ai  shown  in  Fig.  l  3.  The  rotor  consists  of  4  flywheel  fitted  to  a  shaft  and^to 
simplify  the  analytical  worl^lt  will  be  assumed  that  the  rotor  mass  is  con¬ 
centrated  at  a  point  situated  half-way  between  the  bearings  and»  furtheri  that 
the  bearings  are  identical.  The  center  of  gravity  of  the  rotor  is  offset  a  dis¬ 
tance  6  from  the  center  of  the  shaft. 

For  the  amplitudes  of  vibration  of  the  center  line  of  the  shaft  at  ths 
points  1,  2  and  3  we  will  use  the  symbols  e^*  02  and  e^i  rsspeetivsly.  The 
actual  bending  deflections  of  the  shaft  will  be  noted  by  y.  The  relation  be¬ 
tween  the  displacements  at  various  points  on  the  shaft  will  be  as  indioated  by 
the  displacement  diagram  of  Fig:ure  14.  The  amplitude  of  whirl  at  the 
point  2  is  partly  due  to  the  displacement  of  the  shaft  in  the  bearings  and  partly 
due  to  the  deflection  of  the  shaft.  Further,  it  will  be  seen  that  the  amplitvide 
of  vibration  of  the  center  of  gravity  of  the  rotor  is  the  vector  sum  of  9^ 

Figure  15  shows  the  force  diagram  for  the  rotor.  Let  us  denote  the 
angular  velocity  of  the  shaft  by  w  aivd  ths  centrifugal  force  acting  on  the  rotor 
will  be  given  by  the  following  es^ressien: 

C.F.  ■  m  w*  (02  +  4)  <23) 

Since  there  is  no  damping  in  the  shaft,  the  phase  angle  of  the  shaft 
deflection  must  be  the  same  as  that  of  the  applied  centrifugal  force  and  the 
vectors  representing  these  two  quantities  must  be  parallel-  As  shown  by  the 
force  diagram,  the  centrifugal  force  if  composed  of  the  two  components 
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m  w'  £  and  m  w*  to  th«  amplitiidea  6  and  ratpactively. 

The  centrifugal  force  acting  on  the  rotor  Is  equally  divided  between  the 
bearlnga  at  1  and  3  and  the  force  acting  on  each  of  these  bearings  will  there- 
(oti  'iti  equal  to  C.  F./2>  Sines  the  bearings  are  assumi^  to  be  Identieali  the 
■haft  displacements  must  be  the  same  In  both  beadnfs.  We  will  therefore 
drop  the  suffices  and  write  e^  e  <*  e. 

Since  the  forces  applied  to  the  bearings  are  of  constant  nugnitude  and 
rotate  with  a  conitant  angular  veloclty»  the  relation  between  force  and  dis¬ 
placement  can  be  determined  from  the  constant  load  characteristics  of  the 
bearings.  From  Equation  (i)  It  la  apparent  that  introducing  a  value  for  o 
makes  no  difference  to  the  principle  of  the  solution  of  Uie  equation.  Con- 
•equentlyi  the  use  of  curves  obtained  from  the  constant  load  calculations  can 
be  extended  to  cover  the  case  of  a  constant  rotating  load. 

If  we  consider  the  disturbing  force  as  the  unbalance  force  then  it  rotatee 
with  the  eame  angular  velocity  as  the  Journal.  1.  e. .  a  ^  w.  Making  use  of 
the  load  number,  f  (Tah}ee  I  and  IT),  and  taking  into  account  the  rotation,  the 
force  C.  F.  will  be  given  by  the  following  expression: 

*>  ^  u  «  L  D  (K/O*  f  e^  (24) 

For  convenience,  the  attitude  angle  ^  has  been  chosen  to  be  positive 
when  the  force  on  the  fluid  film  leads  the  displacement,  i  -  e. ,  for  negative 
load  nxxmbere.  Here,  f  u/',^  y  are  a  function  of  the  eccentricity  ratio  e  as 
given  in  Tables  1  and  II. 
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L«t  ui  introduce  the  eymbol  ^  for  the  phaee  angle  difference  between 
veetori  fi  and  4  6)  ae  ehown  in  Figure  14. 

By  reeolring  the  dieplaeement  veetora  in  the  x  and  y  direotlon.  we 


arrive  at  the  following  two  equatione: 

e  tin  ^  t>  £  eln  B 


6  4  e^  ^  e  eoa  ^  4  fi  eoe  3  4  y  { 

The  deflection  of  the  shaft  ia  obtained  from  simple  beam  theory  and 


{25a) 

{25b) 


la  given  by: 


c.  F.  l» 

48  El 


Making  use  of  Equations  (23;.  (24)  and  (26),  Equation  (25b}  can  now  be  written; 


)  s  e  eoa  ^46  eoa  B 


Tb4  unknown  3  ean  be  eliminated  between  Equations  (25a)  and  (27)  and 
as  a  result  we  obtain  the  following  relation  between  the  frequency  and  the 
amplitude  of  whirl: 


For  a  further  study  it  will  be  convenient  to  introduce  in  the  above 
equation  a  new  variable  defined  by: 

^  ^  -  «rir 


This  leads  to  a  simple  quadratic  equation,  the  solution  of  which  is: 


1/2/ 


^  ‘  fir^(l)  [•*  - **[)  ‘ 
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M&klag  uie  of  Table  I  axul  the  above  eq^tldai  valuee  can  be  found  for 
7  ai  a  function  of  any  eccentricity  e. 

The  frequency  w  ie  defied  in  terme  of  7  by  Equation  (29)-  Rearranging 
and  eolving  for  u  '7e  get  the  following  formula: 

24  El  7  ^ 

hi  B  -  W  + 


vm  FLPID  riLM  AND  RQtOR  RESONANCE 

Fluid  film  raionanee  it  a  costtqutnct  oi  tiit  quati-eUtt^c  pr<9trti«t 
of  Um  oil  film.  The  phenomenon  may  therefore  be  itudled  by  ae turning  that 
the  rotor  le  rigid.  The  conditiona  of  a  rigid  xntor  may  be  timulattd  by  tatting 
I* /SI  B  0  in  Equation  (29).  Thua*  we  obtain  the  foUowlng  relation: 

7  a  l/nua 

By  tubatltuting  the  exprettion  for  7  in  the  above  equation  we  obtain 
the  following  equation  for  w: 


By  eonaldering  the  foreea  acting  on  the  Journal  it  la  poaaible  to  derive 
the  frequency  reeponae  eharacterietiea  £6r  the  ayatem  and  predict  the  variation 
in  as^litude  of  *fluid  raaonance^  with  apeed.  Referring  to  Fig.  15,  we  ahall 

fcrco  by  th*  ymt^r  P.  the  inetentaaeoue  ooeition  of  the 
Journal  by  the  vector  ti  and  the  force  in  which  the  fluid  film  exerta  on  the 
Journal  by  F.  The  inertia  force  acting  on  the  Journal  will  be  given  by  me* 

Zf  the  Journal  la  vertical  and  there  are  no  ether  external  forcee  acting 
we  can  osqireea  the  conditiona  of  force  equilibrium  for  the  Journal  by  the 
following  vector  equation. 

me  +  F  »  P  (31) 


Let  m£  denote  the  out-of«balanee  of  the  rotor  and  w  ita  angular  velocity- 
The  dleturbing  force  acting  on  the  rotor  will  then  be  given  by: 

_  .  j  i  O  +  ut) 

Pa  ra  0  w  e 


(32) 
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Where  0  ii  the  phaee  Angle  of  the  disturbing  force  relative  to  the 
displacement  vector  e. 

If  we  assume  that  the  journal  center  moves  In  a  circular  path  around  the 

bearing  centeri  we  can  write 

i  ut 

e  °  e  e 

o 

where  e  is  constant  over  any  complete  revolution. 

0 

The  Inertia  force  acting  on  the  Journal  can«  thareforei  be  expressed  by 


.«  a  fut 

me  B  .mu*e  e 
o 


(33) 


If  we  assume  that  a  single  bearing  Is  taking  the  total  load,  Eqxiation 


(30)  bec  -mew 


r  a  H  w  XJ?  (R/O*  f  e* 


(34) 


Making  use  of  Egoations  {32)»  (33).a'  d  (34)  in  Equation  (35),  we  obtaia: 

a  iut.  1  i(uC4  3) 

-m  w*  e^  e  +  M.w  LD  (^)  f  e  ■  m  d  w*  e 

The  real  eompensnts  in  the  above  equation  give: 

•m  w  e  4  iiu  1^.(R/C)*  f  cos  ■  m  6  u  cos  0 

O  6V  ' 

and  from  the  imaginary  components  we  get: 

~  U,  LD  (R/C)*  f  sin  ^  s  m  6  w  sin  0 

Eliminating  0  between  the  above  two  equations  and  rearranging  we 

finally  obtain  the  following  equation  in  w: 

m  U  LD  (R/C)*  f  cos  jT)  [iiLD  (R/C)*  f  J 

—  }  w  +  - 


w“^m  e^)*  -  (m  5)* 


I  WM  I  I 

4v*  (me  )* 

L  o 


(mfi)* 


rf 


a  0 


condition  for  retonnneo  !■  Uint  th«  two  root!  of  tlio  nbovo  aquation 

coincide.  Thle  will  lead  to  the  following  reUtiont 

cos*  ^  n  1  -  ( 2^  )  (35} 

0  ■ 

and  the  reeonance  frequency  will  be  given  by 

m  e  U)  (R/C)*  f  cos  4 
o  ' 


It) 


(36) 
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IX  EQUATIONS  OF  SMALL  OSCILLATIONS 

In  Sections  III  and  IV,  tha  fluid  film  forces  that  act  on  the  Journal  were 
obtained-  We  will  now  consider  the  equations  of  motion  for  small  oscillations 
about  a  position  of  equilibrium  corresponding  to  sonve  particular  external,  steady 
load,  making  use  of  these  forces- 

Introduce  fixed  rectangular  (x.  y)  axes  as  in  Fig.  1 6.  Suppose  now  that 
jinder  a  certain  external  steady  load  on  the  Journal,  it  is  running  at  eccentricity 

ratio  which  corresponds  to  position  A  ,  which  we  choose  as  the  direction 
of  positive  y.  Then  Equation  6  yiwids  for  the  steady  lead  on  the  Journal 


r,  •  0) 

r  r  o 

“  ^t  <‘o' 


(37) 


Consider  next  a  snail  (variable)  deflection  da*  da  of  the  Journal  center  to  A, 
near  the  position  A^*  and  1st  Xf  Y  be  the  added  forces  on  the  Journal  in  the 
X.  y  direction  resulting  from  this  displacement.  Erom  Equations  8  due  to 


changes  in  magnitude  of  F^.  we  obtain  (sinee  u  is  constant) 


df  df 

dF^  ■  y  2  dh  f^  >  u  de  +  ds*) 

*  h£  df 

dF^  ■  \  mZiA  f^  *  w  ds  +^jr 

^  « 

^  df  df  bf 

where  f^.  ^  are  evaluated  at  B  «  e^,  e*  ■  0. 


These  forces  contribute  directly  to  Y.  X  respectively.  In  addition, 
there  are  also  contribiUions  to  X,  Y  due  to  the  directional  changes  da,  given 
by 

Fda  «  -Xufda 
r  r 

-F^  de  *  -  X  u  f^  da 


-  (39) 
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wher«  f  ,  f  are  alto  evaluated  at  E  •  0.  Hence 
r  t  o 

“t  .  1 

SF- 

■'  (40) 

H  b! 

Y  8  X  -  u  da  ^  2dd  f ^  -  u  (^  — ^‘de  + 


uda  f 


2dd  ^ 


Since  £  .  f  are  poeitive  and  Increating  funetlone  of  e  at  well  at  of  s',  both 

'  n  n 

and  tha  vmrleua  darivattvaa  are  poaltlve. 

Let  the  eomponentt  of  the  dltplaeement  A^A  along  the  Xi  y  direction 
be  (>  T]  (tee  Fig.  16  )  and  let  a  be  meatured  from  y-axlt.  Then  for  email 

4,  T)  and  a,  there  results  from  Fig..  J6. 


de  o 


da  B 


t  C  ’ 
o 


de'  ■  —  *» 

u>  Cw 


"  ■  ri" 


Equatlont  (40)  yield  , 

„  .  .  r  “'r< 


X  ■  X  . 


df  • 


Y  B  X 


r«r5  zt  i  bt  ^ 

t  ,  r  rfl  r’l 

•r^  *  ITT  ■  “‘sr  c  ■  sF’c 

Loo  j 


Because  of  the  assumption  of  small  a,  Equations  41  and  42  are  not  valid  for 
small  £ 

Now,  suppose,  as  in  Fig.  17,  that  a  rotor,  consisting  of  a  single  disk  is 
symmetrically  mounted  on  a  shaft,  which  is  carried  in  two  similar  and 


symmetrically  placed  bearings.  Let  k  be  the  shaft  stiffness  constant  We 


assume  that  the  deflections  of  the  shaft  in  each  bearing,  both  for  the  steady 
load  deflection  and  for  small  further  deflections  are  symmetrical  about 
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the  middle  plaune  of  the  rotor,  end  that  eteady  load  and  the  added  loreee 

X,  Y  on  each  Journal  are  equal .  The  ahaft  eegment  to  the  left  of  the  rotor 

hai  the  force  X,  Y  4Y  exerted  on  it  by  the  lubricant  on  ite  Journal,  and  muet 

0 

have  an  equal  and  oppoelte  force  exerted  on  It  by  the  rotor  mail,  (which  alio 
exert!  a  elmllar  force  on  the  ihaft  eegment  to  the  right)*  Thui  the  ihaft 
center  deflect!  relative  to  the  Journal  center!  an  amount 

X  «  -  2X/k,  y  ■  -2Y/k  (43) 

giving  rlee  to  a  net  dliplacement  of  the  rotor  center 

X  +  ^,  y  +  ^ 

We  euppoce  that  the  eteady  force  •>2Y^  le  an  external  load  on  the 
rotor  (gravity).  The  remaining  forcee.  2X  and  2Y.  aaauming  that  the  rotor 
l!  balanced,  can  only  be  uied  ae  *inertla  forcee^.  Hence 
2X  •  m  (if  +  (S)  >  .  kx 

1451 

2Y  ■  m  (9  +  'i)  ■  -  ky 


and  eubitltutlng  from  Equation!  (42)  we  get 


wf  5  2f,  e  A 

'  -T§--TV*-sr^i*3r3 


f  (i%<)  • 


wf-t  2f  i  wdf 
t  ,  r  r 


>  .  — -y.  a  — 

e  C  e  C 
L  o  0 


_ r  n  _ r  n 


^(y  +  *ti)  H 


Here  the  function!  f^,  f^  and  their  derivative!  with  reepect  to  e,  e* 
are  all  evaluated  at  the  equilibrium  eccentricity  ratio,  and  lor  e*  «  0.  They  are 
given  in  Fig.  18 '21  for  the  ineomprefilble  eaaee  and  in  Fig.  22-25  for  the 
compreiilble  caeae. 
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The  dli£ere&ti«l  eqiutlont  (46)  are  linear  in  the  variable  (.  t);  xi  y  and 
■olutione  of  them  esdet  in  which  time  entera  aa  an  exponential.  Theae  may  be 
expreiiad  in  dlmenalonleae  form  at 


VT 
e  , 


where  T  ■  u  t, 
o 


(47) 


Hera  ia  the  critical  apeed  of  the  rotor  mate  m  with  a  ahaft  atiffnaaa  k. 
From  the  right  hand  equation  pair  of  equation  (46  )•  there  raaulta 


X  ®  - 


m  (vw^)* 
k  +  m  (vtti 


m  (V w  )* 

^  ’  k  4  m  (V  w  )*  ^ 

o 


Oniet  of  IniUibillty 

It  m&y  be  ihown  th&t  for  imAll  valuei  of  Bi  all  sljt  roots  of  Equation  (51) 

He  In  the  lefrbalf  of  the  complex  v  plane  This  means  that  the  solutions  of  the 

small  displacement  equations  (46).  about  the  position  of  equilibrium  coTre8<^ 

ponding  to  the  given  constant  load  on  the  rotor.,  are  stable.  On  the  other  hand. 

if  s  is  large,  it  can  be  shov  n  that  some  of  the  roots  He  in  the  real  half  of  the 

complex  plane,  so  that  the  ttolutlons  of  Equation  (46)  become  unstable.  The 

transition  from  stability  to  inNtablllty,  one  may  further  show,  can  never  take 

place  by  a  root  crossing  at  v  =0.  Therefore,  this  transition  can  only  occur 

as  a  result  of  roots  crossing  the  pure  Imaginary  v  axist  as  illustrated  in 

Figure  26.  We  t^w  proceed  to  find  when  this  crossing  takes  place.  (Note 

first  that  to  study  the  roots  of  Equation  (51),  it  is  convenient  to  keep  constant 

the  value  of  t  at  which  f  and  T  and  their  derivatives  are  evaluated. ) 
o  r  t 

Writing 


,  .  E  V* 

'  ■  rr^ 

Equation  (51)  becomes 


sf^+2vf^+eC 


(5  3) 


df. 


-8  f  +  2  V  f 
t  r 


**5r 


Sf  bt 

“■sf  ' 


(54) 


To  find  the  threshold  of  instability,  we  seek  pure  imaginary  roots  of 


(51).  For  such  a  root,  both  the  real  part  and  the  pure  imaginary  part  of  f  (v) 
in  (51)  must  vanish  separately. 
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EquAtlng  to  sero  the  imagin&ry  part  o£  Equation  (54),  we  have; 


it 

1 

c« 

*'i  •'sr 

•  'r  ♦  ‘  '  S?' 

+  V 

b{ 

b{ 

-S  f  V 

t  dE* 

Since  V  ^  0: 

.  r  df  ]  Sf  bf 

+  **•<*, -sf  IT  ’  '  ° 

This  is  linear  In  (*  and  yields  for  ^/s  an  explicit  function  of  t  which  turns  out 


negative: 


s 


bt  bf  bt  b{ 

- ^ - g- - 


Next,  equating  to  aero  the  real  part  of  Equation  (54),  we  have: 


dfj 

5f 

•  ^  ^  ^  •IT 

f  --L 

t  Se* 

b{^ 

+  2  V* 

bf 

‘r 

Upon  dividing  by  s*.  substituting  for  C/s  from  (55),  and  solving  for  (v/g)*, 


there  results: 


This  Is  a  function  of  e  only,  and  turns  out  to  be  positive,  yielding  a  pure 
imaginary  value  of 

V  =  1  v‘ 


where 


(57) 


a 


±\fU^ 

Ct 


<8  the  ratio  of  the  oacillatlon  frequency  of  the  journal 


The  quantity 

center  to  the  shaft  running  frequency  at  the  limit  of  stability.  Note  that  it  is 
independent  of  E  and  depends  only  on  e. 

We  now  return  to  Equation  (5  3'  and  put  it  in  the  form: 

5 


\s\ 

V 

<  -  B  E 

V 

• 

(58) 


Recalling  (55),  (56)  and  dividing  (58)  by  f ^  (e)  '  obtain: 

8  E  1  .  „ 


8  + 


fj  (€) 


f2(e) 


(59) 


The  solution  a  (-  u/ij^)  of  this  quadratic  equation  depends  both  on  £ 
and  on  £  and  is  thus  dependent  (see  Equation  52)  on  the  bearing  geometry, 
fluid  viscosity,  rotor  mass  and  shaft  stiffness.  The  value  of  s^obtained  by 
solving  Equation  (59)^expres8es  the  rotor  speed  at  the  limit  of  stability  as  a 
function  of  the  masS’Shaft  critical  speed. 

It  should  be  kept  in  mind,  however,  that  the  critical  speed  of  a  mass- 
shaft-bearings  system  should  be  calculated  using  the  combined  stifiness  '>{  the 
shaft  and  lubricant  film.  For  a  symmetrical  2>bearings  system,  this  critical 


speed  may  be  calculated  as  follows: 


a)  Sh&ft  Stiffness  =  k 

b)  Lubricant  film  stiffness  =  H-  ^  ^ 

de  V  C  de 

sk  df 
*  2E  d£ 


The  critical  speed  of  the  system  is  then: 


u 


CR 


m  ( 


1  +  i  . 

k  2  ks  /  de 


Its  ratio  to  the  shaft-mass  criticali  u  ,  is 

o 


In  defining  the  lubricant  film  stiffness  as  above,  it  is  customary  to  use 

the  radial  film  stiffness  dF  /de. 

r 

Now,  the  ratio  of  the  rotor  speed  at  the  limit  of  stability  to  the  critical 
speed  of  the  rotor- shaft -bearings  system  can  be  obtained* 


The  value  of  s  used  in  this  equation  being  that  obtained  from  the  solution  of 
Equation  (59). 
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Example: 

Coniider  the  rotor-ehAft-beAringa  ayatem  of  Figure  27. 

R  B  1  inch 

^  =  0.5  And  I 

£  .  10- 

k  a  1.25  X  10*  lb. /in.  (From  Fig.  27  for  a  ateel  rotor  4  ahAlt) 

m  a  0. 1  lb.  -  tec.  (From  Fig.  27  for  a  ateel  rotorakehAft) 

The  dimenalonleaa  number  E  (defined  in  EquAtlon  52)i  la  a  function  of 
bearing  geometry,  ehaft  atiifneea  and  fluid  vlteoalty.  In  the  example  worked 
out  here,  the  rangetO.  1^  £d.l00  la  covered.  For  the  rotor  ahown  In  Fig.  27 
thia  eorreaponda  to  a  viacoaity  ranging  from  order  10*^  to  order  lO’^ 
lb.  -ICC.  */ln.  The  range  of  eeeantrlclty  ratioa  covered  in  thia  eaEample  la 
O.ldeSO.e, 

The  film  atiffneaa  uaad  In  calculating  »he  natural  frequency  (per 

Equation  34)  la  that  due  to  the  radial  reatdring  force.  The  natural  frequenelea 

calculated  below  will  then  be  the ae  along  tha  radial  direction  (u  )  . 

Sf  af  Sf  af 

The  value  a  of  f^.  f^i  *5“'  obtained  from  the 

aOlutlOn  of  the  dimenalonleaa  Reynolda  Equation.  Theae  are  given  in  Table  I 
and  Figurea  18  through  21  for  the  incomp  re  aaible  caaea  and  Table  U  and 
Fi^.:rea  22  through  25  for  the  compreaaible  caaea. 

Introducing  theae  valuea  into  Equationa'§2,  55,  59.  and  60.  we  obtain: 


(for  the  incomp  re  aaible  caaea); 
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Bearlng  L/D  =0.5 


e 

E 

1 

a 

1^1' 

0.  2 

0.  1 

-  0  9429 

-0. 1411 

0.5 

0.  1 

-  3. 1968 

-0.  1252 

0.8 

0.1 

-13.4452 

-O.OSS43 

0.2 

1.0 

-  0.9429 

-0. 1411 

0.  S 

1.0 

-  3.  1968 

-0. 1252 

0.8 

1.0 

-I  3.4452 

-0.05543 

0.2 

10. 0 

-  0.9429 

-0. 1411 

0.5 

10.0 

-  3.  1968 

-0.  1252 

0.8 

10. 0 

-13.4452 

-0.05543 

0.2 

100.0 

-  0.9429 

-0. 1411 

0.5 

100.0 

-  3.  1968 

•0. 1252 

0.8 

100.0. 

•  13.4452 

-0.05543 

I  jv\  u 


0.  3756 

2.6096 

0.9852 

2.6488 

0.  3538 

2.8104 

0.9974 

2.8177 

0.2354 

4.2438 

0.9998 

4.2446 

0.  3756 

2.  1842 

0.8533 

2.5597 

0.  3538 

2.67^3(9 

0.9741 

2.7450 

0.2354 

4.2105 

0.9985 

4.2168 

0. 3756 

0.6307 

0.2723 

2.  3162 

0.  35  38 

1.6659 

0.7325 

2.2743 

0.2354 

3.8919 

0.9839 

3.9556 

0. 3756 

0.06678 

0.02910 

2.2948 

0.  3538 

0.2533 

0. 1 315 

1.9262 

0.2354 

1.9267 

0.7751 

2.4857 
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Bearing  L/D  »  1 


e 

E 

1 

s 

(■f 

>1  ■ 

s 

^“CR^  r 
*^0 

w 

^“cR^r 

0.  2 

0.  1 

•  2.6761 

-0. 1377 

0.37 

2.6763 

0.9948 

2.6903 

0.  3 

0.  1 

.  4.324 

•0, 1338 

0.37 

2. 7224 

0.9971 

2.7303 

0.7 

0.  1 

.18.484 

-0.0846 

0.29 

3.4179 

0. 9998 

3.4186 

C.  8 

0.  1 

.29.2005 

-0.04898 

0.22 

4.5169 

0.9999 

4.5173 

0.2 

1.0 

-  2.6761 

-0.1377 

0.37 

2. 5146 

0.9488 

2. 6503 

0.  3 

1.0 

-  4.324 

-0. 1338 

0.37 

2.620 

0.9704 

2.7000 

0.7 

1.0 

-18.483 

-0.0846 

0.29 

3.3936 

0.9974 

3,4024 

0.8 

1.0 

•29.2003 

.0.04896 

0.22 

4.5016 

0.9992 

4.5052 

0,2 

10. 0 

-  2.6761 

-0.1377 

0.37 

1.4108 

0.5798 

2.4332 

0,  3 

10.0 

-4.324 

-0.1338 

0.37 

1.8120 

0.7263 

2.4950 

0.  7 

10.0 

-18.483 

-0.0846 

0.29 

3.1603 

0,9737 

3.2457 

0.8 

10.0 

-29.2005 

-0.04898 

0.22 

4.  3506 

0.9917 

4.  3872 

0.2 

100.0 

•  2.6761 

-0. 1377 

0.  37 

0. 1951 

0.0834 

2.  3420 

0.  3 

100.0 

-  4. 324 

-0.1338 

0.37 

0.3198 

0. 1390 

2.300 

0,7 

100.0 

-18.483 

-0.0846 

0.29 

1.6553 

0. 7000 

2.3647 

O.B 

100.0 

-29.2005 

-0.04898 

0.22 

3.1199 

0.8998 

3.4675 

The  ratios  .  ■  'r  '  for  these  (incompressible)  cases  are  clotted  in  Figure  28 

rad 

as  a  function  of  the  oeeentrieity  ratio,  L/D  and  dimensionless  number  E.  The 
figure  indicates  ioac  wiitie  for  iow  eccentricity  ratios  the  instability  occurs  at 
approximately  twice  the  critical  speed,  this  number  increases  with  increase  in 
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eccentricity  ratio.  Thus  the  onset  of  Instability  for  eccentricity  ratios  of 
0.  8  is  around  three  to  four  times  the  critical  speed.  This  conclusion  agrees 
with  observations  which  show  that  stability  Increases  with  increase  in 
eccentricity  ratio.  However,  it  also  shows  that  Instability  may  occur  even  at 
very  high  eccentricity  ratios.  These  have  not  been  observed  for  very  few 
machines  in  practice  operate  at  speeds  much  above  twice  the  first  critical. 

The  number  —  |“| ,  shown  in  the  above  table  (where  i  *  V^)i  represents 
the  ratio  of  the  frequency  of  the  oscillation  of  the  shaft  center  to  the  running 
frequency  of  the  shaft,  calculated  at  the  onset  of  instability,  as  discussed  on 
page  35.  This  ratio,  which  was  shown  to  be  independent  of  E,  is  plotted  in 
Figure  29  as  a  function  of  the  eccentricity  ratio.  Note  from  this  figure  that 
for  a  wide  band  ot  eccentricity  ratios  (between  0. 2  and  0.  5),  the  frequency  of 
the  shaft  oscillations  at  the  onset  of  instability  is  between  35^  and  40^  of  the 
running  frequency  of  the  shaft. 

Calculations  were  also  made,  using  the  analysis  described  above,  for 

bearings  with  an  L/D  of  1  and  using  compressible  theory.  The  forces  and 

gradients  used  in  these  calculations  are  those  given  in  Table  II  and  Figures 

22-25.  From  these  calculations,  the  ratio  - - -  is  plotted  in  Figure  30 

rad 

as  a  function  of  eccentricity  ratio  and  for  values  of  E  -  10  and  iOO.  Note, 
from  this  figure,  that  for  the  compressible  case  loo,  the  limit  of  stability 
increases  with  increase  in  eccentricity  ratio  At  low  eccentricities,  however, 
instability  sets  in  at  very  low  speed-  In  fact,  Figure  10  suggests  that  at 
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eero  eccentricity  ratio,  oscillations  will  be  unstable  at  all  speeds.  Some  ex¬ 
perimental  verification  of  this  is  given  in  Reference  4,  where  it  la  shown  that 
for  vertical  rotors  (e  a  o)  instability  exists  at  all  speeds. 

It  must  be  pointed  out  here  that  the  mesh  sice  used  {in  the  digital  com¬ 
puter  calculations  of  the  fluid  film  forces  and  gradients)  for  the  incom¬ 
pressible  Cases,  was  approximately  three  times  as  coarse  as  that  used  for 
the  compressible  cases.  The  reason  for  this  is  the  much  greater  complexity 
of  the  non-linear,  compressible  fluid  Reynolds  equation.  Checks  were 
made  in  a  few  calculations,  using  a  finer  circumferential  grid  and  these 
showed  that  the  values  of  total  dimensionless  force  (f  obtained 

with  good  accuracy  in  both  the  compreesible  and  incompressible  cases. 

j  f 

However,  the  calculated  values  of  the  attitiide  angle  s  Tan*  )  may  be 

*r 

in  error  by  two  or  three  degrees  in  some  of  the  compressible  cases. 

Since  the  relative  values  of  the  tangential  and  radial  components  of  fluid  film 
forces  are  very  critical  In  stability  investigations,  it  is  planned  to  review 
the  results  given  in  Figure  30  using  more  accurate  computations. 

The  difficulty  in  obtaining  accurate  numerical  solutions  of  the  non¬ 
linear.  compressible  fluid  Reynolds  equation  Illustrates  the  advisability  of 
evaluating  the  order  of  magnitude  of  each  term  of  the  equation.  Thle  allows 
those  terms  whose  order  of  magnitude  can  be  shown  to  be  small  (relative 
to  the  other  terms  of  the  equation),  to  be  omitted  from  the  equation  before 
numerical  solutions  are  obtained  for  a  wide  range  of  cases.  Failure  to  do 
this  greatly  increases  the  difficulty  and  computation  expense  required  to 


obtain  accurate  results. 
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X  EQUATIONS  OF  MOTION  TOR  LAAOE  DISPLACEMENTS 

We  now  eoniider  the  general  equetione  of  motion  of  the  Journal  center 
and  the  rotor  center,  without  making  the  aeeumption  uied  in  lection  IX  that 
thie  motion  la  cloee  to  the  ateady-etate  poaition. 

Suppoae  that  under  a  certain  applied  ateady  load  on  the  rotor,  euch  aa 
gravity,  the  Journal  center  la  deflected  to  Q^,  corraapondlng  to  e  ■  e^,  and  that  co¬ 
ordinate  axea  are  choaan  aa  in  Figure  11  lo  that  y  axia  paaaea  through 


Then  the  coordinate  a  of  are 


^  B  0 


n  “  C 


(61) 


The  foreea  exerted  by  the  lubricant  on  the  Journal  at  each  bearing  (half  the 
applied  ateady  load)  are  given  by 


F  ^  o  X  «  f.  (e^,  0), 
xo  to 


F  »  -X  wf  (e  .  0) 
yo  r  o 


(62) 


Theae  force!  alao  cauae  a  ateady  deflection  Q^A^  of  the  abaft  center 


xjid  rotor  center  of  amount! 


o  xo 


y  ■  -2  F  /k 
’^o  yo 


(63) 


Conaider  now  a  general  motion  of  the  Journal  center  to  the  point  Q 

apecified  by  the  polar  coordinate!  (eC,  a).  The  Journal  center  ia  now  at 

5  ■  C  e  aln  e,  i)  ■  C  e  coa  o  (64) 

The  force!  on  the  Journal  in  the  radial  and  tangential  direction!  are  (from  Eq.  8 
and  8a) 

E  B  m\  -  2o)  f  (c,  C  ) 

r  '  (65) 

F^  a  X  («  -  2a)  f^  (E,  e’) 


where 


e'  B 


(1  -  2a') 


(65a) 
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and  along  x.  y  diractiona 

m 

■  *  X  (w  -  2a)  (e,  t»)  tin  »  +  (Si  **)  cot  a 

f  ' 

*'y  ■  •>(«*■  24)  {e,  V)  cot  tt  +  (*•  *')  •!»  a 

»  , 

The  ihalt  center  deflection  it 


2X 

*  "  ■^V'  ■  (w  •  Zk)  -f^  tin  a  +  fj  cot  a 

•  « 

2x  r 

y  B  -2F^k  a  (w  -  24)  f^  cot  a  +  f^  tin  a 

L  * 


and  the  net  deileetiona  of  the  rotor  center  are: 

2X 

4  +  X  a  C  e  tin  a  +  (w  •  24)  tin  a  •  f^  cot  a 


n  +  y  ■  C  e  cot  o  +  ^  (w  -  24)  f  cot  a  *  f^ 


tin  aj 


Subtracting  the  ateady  force  componentt  Equation  (62)  from  Equation  (66) 


ere  obtain: 


■  -X  u)  fj  (e^i  0)  ♦  X  (u  •  24}  -f^  tin  a  ♦  f^  cot  a 

r  -  m 


r  -  F 

y  yo 


r  a 

(e  f  0)  •  X  (v  *  24)  f  cot  a  ♦  f  tin  a 

r  w  T  t 

'  ■ 


Thete  foreea  con  only  be  uttd  up  in  producing  acceleration  of  the  rotor  matt. 


That. 


T  '*  ♦  *!’>  "  T  'i'  ♦  ''i>  *  *■»•»■. 


or  more  explictly, 

[  ^  jjir  +  1  j  (w  -  2o)(-f  tin  a  4>  f  cot  a) 


m  d' 


y  yo 


"  <,  (e,>  0)  ♦  it  ain  a) 


FS*’*  'll  ‘r”"*  I  “  u  f,  (e^.  0)  .  ^  |L(e  cPi  a) 
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Maitli 


Rotor 


If  w«  eontider  tht  cate  of  mataloiit  rotofi  tht  taraa  involvlai  m  can 
ba  dalatad  from  Equation  71  which  than  raducaa  to: 


(■-¥ 

* 

- 

-f^  (E. 

s')  ain  E  4  f^  (e.  I'}  coa  a 

« 

(■•a 

f  (e, 

7 

■ 

e*)  coa  a  4  f^  (e»  i*)  ain  a 

\ 

. 

■  f  {e 0) 

I  0 


a  f  (B  ,  0) 
T  O 


(72) 


Theaa  are  differential  equationa  from  which  d.  i  can  ba  aolved  aa 
functiona  of  e»  a. 

The  ateady  atate  aolutlon  ia  given  by  s  a  a  a  o.  Equationa  72 
may  be  tranaformed  by  multiplying  firat  by  -  ain  Bi  coa  a  and  adding^ 
then  by  coa  a,  ain  a«  and  adding.  There  reaulta 


24 

1  - 

w 


£,  (e,  ?•)  a  -  f  (e  ,  0)  ain  a  4  f  (e  .  0)  coa  a 
f  to  r  o 


f.  (e.  e')  a  r  (e^t  0)  coa  a  +  f  (c^,  0)  ain  a 

X  to  TO 


(73) 


Introducing  aa  in  Figure  32,  the  reaultant  force  f  f^ven  in  dimenaionleaa 


form 


f  (e,  e')  a  ^fj  +  f J  (74) 

and  the  direction  angle  ^  that  thia  force  makea  with  the  (inner)  radial 
diret 


Kf  a  tan'^  f^/f^ 


(75) 


We  note  alao  the  relation 


f^  a  f  coa  %f, 


f^  a  f  ain  f 


(76) 
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D«nete  by  £^i  the  veluee  of  i  end  i/  et  the  peeltione  of  equilibrium. 

0  o  o  o 

Thui,  by  taking  the  equere  root  of  the  eum  of  iquaree  of  Equetiohi  1)  end 
their  ratio,  one  o^  -  the  following  equivalent  form  for  theee  equatlone; 

f  ie.  !•)  »  f^  (78) 

tan  %/  B  tan  (a  +  (79) 


Hence 


“  o  +  (80). 

Introduce  the  parameter  t  by  meane  of 

T  B  Wt  (81) 

Then  Equation  (78)  can  be  replaced  by 

0  -  2a‘)  f  (e.  i')  «  f  (e  ,  0)  b  f  (82) 

o  o 

By  meane  of  plote  of  f  and  i/  in  the  e  -  plane,  euch  ae  thoee  of 

Figuree  5  and  6 ,  the  locue  of  a  ehaft  center  initially  dleplaeed  from  Ite 

equilibrium  poeitlon  (e  .  !('  )  to  eome  other  position  (e,  i/)t  can  be  obtained 

0  o 

numerically. 

The  equations  required  for  thle  and  a  description  of  the  procedure 
follow: 

From  Equation  (80)  we  have 

a  B  ^  ^  (8  3) 

o 

and  from  Equation  (82)  we  have 


(84) 


For  a  •mall  dimenslonleta  time  increment 


a  "I  J.1  -  T 

n+l  n 


we  can  write 


•  E  +  e*  •  iji 

T*T  TBT  TBT 

n+l  n  n 


L  o  a  +  a*  •  ^  T 

■t»T  TBT  TBT 

n+l  n  n 


Now.  coneider  a  shaft  center  whose  static  equilibrium  position  is 
(e  .  0)  and  let  it  be  displaced  to  a  position  (e  a  )  at  the  instant 

O  TBT  >  TBT 

O  O 

T  ■  T^*  Its  locus  can  then  be  traced  as  followe  (See  Figure  13a)  t 
1.  From  Figures  S  and  6,  obtain  the  dimensionless  force  f  and  the 
angle  which  correspond  to  the  position  of  static  equilibrium 

(e  ■  e*  ■  0). 

0 


From  Equation  8  3,  calculate 


8  p  +  a 

T*T  O  T*T 

O  0 


3.  Locate  the  point  Figure  6  and  hence  read  off  e*  ^ 


4.  Locate  the  point  (e  ,  S*  ^  )  in  Figure  8  and  read  off  f 

TBT  TBT  *  T“T 

0  0  0 

5.  By  means  of  E^quation  (82),  calculate 

Ilf 
a*  o  1  - 

TBT  2  2  f 

o  o 

6.  For  a  small  time  increments  t  a  t^>t^,  calculate  from  Equation  (86) 

a*  B  a  +  a*  •  ^t 

TSTi  TBT  TBT  " 

i  0  0 

7.  From  Equation  (64)  calculate 
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8.  From  Equation  {65a}  obtain 


C*  a 

TOT 

0 


[e*  (I  -  2a«) 

I  Tax 


9‘  Tvom  Equation  (8  3)  calculata 

£  8, .  e  +6*  •  ^  T 

TBT,  XOT^  XBT 

10  0 

Tha  now  location  of  the  abaft  canter  (e  .a  )'ia  now  defined. 

X*Tj  V“Tj 

The  procedure  la  repeated  for  a  new  time  interrmli^  t  "  ^2  -  and  the 
location  of  the  ahaft  center  at  time  t  a  la  almllarly  obtained.-  Stepwiae 
repetition  of  thla  procedure  will  than  define  the  complete  locua  of  the  ahaft 
center.  Figurea  33  (b).  (c)  and  (d)  ahow  three  poaaible  loci. 

In  Figure  33  (b),  the  ahaft  center  leada  barb  to  the  poaitlon  of  atatic 
equilibrium  ao  that  the  condition  ia  atable. 

In  Figure  33  (c).  the  ahaft  center  deacribea  a  cloaed  loop.  Thla  la 
a  condition  of  atable  oacillation.  The  period  of  theae  oacillationc  ia 

In  Figure  33  (d),  the  locua  of  the  ahaft  center  apirala  out  and  the 
oaciUatlona  are  unatable. 

As  an  example  uf  this  procedure,  a  numerical  integration  was  carried  out 
for  the  case  where  the  Journal  center  was  displaced  radially  from  Its  steady 
state  position  at  £  =  0.  l  to  an  Initial  position  at  e  =  0.4.  The  fluid  film  forces 
used  were  those  obtained  from  incompressible  theory  for  a  bearing  with  an 
L/D  of  1 .  The  results  of  this  integration,  carried  through  for  the  first  loop, 
are  shown  in  Figure  34  These  indicate  that,  for  the  conditions  assumed,  the 
Journal  center  spirals  in  towards  the  steady  state  position. 
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XI  EXPERIMENTAL  VEimiCATION 

Scvaral  ebianratioii«  may  b«  ^ted  which  varliy  tha  thaoryi  and  to 
aoma  axtant  clarify  thia  vary  difficult  and  controvaraial  aubjaet. 

Eifura  35  chowa  tha  amplituda  and  fraqoaney  obtaiaad  axparimantally 
on  a  vartieal  700  HP  eannad  rotor  motor.  Tha'firat  ayatam  critical  for  thia 
motor  waa  alao  maaaurad  by  axeitation  and  waa  found  to  ba  at  28. 6-29 
cyelea/aac.  Similar  raaulta  wara  obtained  on  four  othar  unita  of  aame  con- 
atruetion.  Thia  iUuatration  la  given  for  it  eitaa  two  forma  of  inatability 
and  ahowa  that  thara  ia  actually  tranaition  batwaan  tham.  Tha  amplitude 
of  the  Half  Frequency  Whirl  ia  general  ia  quite  email  andi  tharaforei  the 
rotor  fepeed  maybe  ineraaaed  without  eauaing  damage  to  tha  baaringa.  On 
tha  othar  hand,  tha  amplitude  of  Raaonaat  Whip  ia  tM>  large  for  aatiafaetory 
operation.  Raf.  4  givea  a  number  of  axamplaa  of  Half  Frequency  Whirl. 

It  alao  poiata  out  that  eompreaaible  fluid  baaringa  are  alao  auaeeptible  to  thia 
form  of  inatability.  Thiii  point  waa  further  verified  more  recently  in  Ref.  5. 

Figure  35  ahowa  that^dep ending  on  tha  apaed  of  operation,  one  may  find 
different  forma  of  inatability.  It  la  not,  thareforaf  aurpriaing  that  conflicta 
in  obaervation  exlat.  Tha  Half  Frequency  Whirl  haa  not  been  obaarved  too 
often  becauae  ita  amplitude  of  vibration  in  general  ia  araall  and  may  not  cauae 
severe  problema  to  the  equipment.  Thia  form  of  inatability  ia  more  common 
for  vertical  rotora  and  lightly  loaded  machinea  aapecially  of  the  eompreaaible 
fluid  type.  The  Resonant  Whip  on  tha  other  hand  which  occurs  in  general  at 
speeds  twice  the  first  system  critical  has  large  amplitude.  It  has  been  shown 
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however  (Tig.  28  &Ad  30)  thit  thfl  occur&nce  of  thii  ttubllity  ii  dependent  on 
eccentricity  retio  and  the  parirseter  E.  Thui.  once  the  Initability  occuri 
it!  frequency  may  increaie  or  decreaie  with  increase  in  ipeed  depending  on 
the  actujil  position  of  Joum^  center  and  the  value  of  E  at  the  specific  speed. 
For  as  speed  increases,  the  eccentricity  and  the  viscosity  decrease  and  they 
both  influence  the  threshold  of  instability.  If  the  damping  increases  vith 
speed,  the  tendency  is  to  decrease  the  Resonant  Whip  Frequency,  the  reverse 
is  also  true. 

As  stated  previously,  observations  have  air  j  been  made  on  the  inception 

of  instability  with  turbine  or  motor  load.  This  ;s  to  be  expected  for.  ii  we 

refer  to  the  torque  on  a  turbine  or  motor,  we  find  that  there  is  an  eccentricity 

between  the  rotor  and  the  stator.  Thus,  uneven  torques  are  present,  see 

Figure  36.  In  the  case  of  a  turbine,  at  the  point  of  closest  approach  the 

steam  or  gas  leakage  path  is  decreased  and  therefore,  the  force  at  that  point 

is  larger  than  that  180°  away.  In  the  case  of  a  motor  or  a  gt.»..,rator,  the 

flux  increases  and  this  has  a  similar  effect.  Thus,  there  is  a  resultant 

l 

force  at  O'  which  is  at  right  angles  to  the  displacement.  It  acts  in  the  same 
way  as  the  tangential  fluid  film  force  and,  therefore,  must  be  Included  in 
the  stability  analysis. 

In  sections  IX  and  X  It  has  been  shown  that  instability'  is  dependent  on 
the  fluid  film  propertlei  within  the  bearing  clearance.  It  is  therefore 
essential  to  evaluate  these  properties  accurately.  These  sections  show  that 
the  inception  of  instablUty  depends  to  a  very  large  extent  on  eccentricity 


ratio  and  lyatam  critical.  A  fancral  ■tatamant  can  ba  mada  that  tha  thraihold 
of  initability  inereaaei  with  incraaaa  in  aeeantriclty  ratio. 

Slnca  tha  atilfnaia  of  the  fluid  filnn  ia  axtramaly  important  in  pradlctinf 
•yetem  criticala  and  alao  fraquaneiaa  of  initability.  tha  point  daiarvae  a 
brief  diieuieioQ.  It  haa  baan  ahown  in  Raf.  4  and  confirmed  further  in  Ref.  5 
that  Equation  2 1  will  predict  accurately  the  thraahold  of  Half  rraquaney 
Whirl.  Thia  equation,  howarer,  contain!  the  radial  fluid  film  atiffnaaa  which 
haa  to  ba  obtained  accurately  either  by  theory  or  from  teata.  In  Table  in 
of  Ref.  5  good  compariaon  ia  giran  between  the  maaaurr.d  and  calculated 
Whirl  rrequaneiaa.  The  correlation  between  tha  two  ia  further  improved  if 
compraaaibla  theory  ia  uaad.  Figure  37  employe  compra a aibla  theory  and 
ahowa  tha  radial  componest  of  force  vi  eeeantrlelty  ratio  for  bearing  B 
diaeuaaed  in  Ref.  5  (Figure  1  and  Table  3).  For  the  load  eorreapondiag  to 
point  (a)  the  radial  atlifneai  for  the  compreaaible  eaee  ie  actually  lower  than 
that  obtained  from  Ref.  5.  Figure  1 .  Thle  will,  therefore,  lower  the 
calculated  thraahold  of  Whirl.  Tha  reverae  ia  true  for  the  two  higher  loada 
b  and  c. 

Figure  37  waa  plotted  to  ahow  another  relevent  point.  Mote  that  the 
radial  component  of  force  la  a  function  of  both  eccentricity  ratio  and  apeed 

Since  inereaeed  load  produce!  both  %ncreaae  In  eccentricity  and 
threahold  of  initability,  the  queatlon  ariace  aa  to  where  doea  one  calculate 
the  radial  itiffneia,  e.g.  ,  at  b  or  b',  c,  c'  or  e".  It  ii  fortunate  that  the 
atiffneia  doea  not  vary  toe  much  between  point  b  and  b'  or  c  and  c"  aa  ahown 
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in  Figure  3*7.  Further,  elnee  etiffneie  entere  to  the  one  heif  power  in 
Equetion  21,  the  error!  in  threehold  of  whirl  will  be  relatively  emell. 

In  the  eiieuletioi;  of  lyitem  eriticelii  the  fluid  film  itiifneei  muit  be 
coniidered.  In  the  uee  of  hydrodynemic  beehngi,  there  ii  aome  queation 
ea  to  which  itiffneaa  ihould  be  uaed.  For,  ea  ahewn  in  aactiona  HI,  IX  and 
X,  there  ere  both  e  rediel  end  e  tengentiel  atiffneaa  end,  in  generel,  the  two 
ere  not  cquel.  In  rlgoroua  enelyaia  one  muat  uaa  the  derivetive  of  force  with 
reapcet  to  diaplecement.  The  increment  of  force  muat  contain  coupling 
terma.  Since  the  atiffneaa  veriea  with  attitude,  there  will  be  e  number  of 
critieela  which  will  cover  e  band  of  frequanciea.  Large  variation  in  redial 
end  tangential  atiffneaa  will  yield  a  wide  frequency  band.  Convaraaly,  in 
the  caae  where  the  atiffneaaea  are  nearly  equal  in  ell  direetlona,  the  critical 
frequency  band  will  be  narrow. 

In  the  caae  of  hydroatatic  baaringa,  the  problem  of  calculation  of 
critical  apeed  ia  very  much  eimpler,  ee  ahown  in  aection  IV.  For  auch 
bearing!,  the  atiffneaa  aa  a  function  of  load  remain!  eonatant  at  leaat  up  to 
eccentricity  ratio  of  0.  S.  Thua  natural  frequency  la  independent  of  external 
force.  Tbla  aimplifiea  the  calculation  of  the  natural  frequency.  Experimental 
verification  of  calculated  natural  frequency  waa  made  recently  on  a  Hellum/Nitro- 
gen  compreaaor  wboae  rotor  maaa  waa  Zx  10*^  lb.  -aee.^/in.  and  whoae  rotor 
atiffneaa  waa  3x10^  lb.  /in.  The  measured  natural  frequency  with  100  pala 
aupply  preaaure  waa  698  cyclea/aec. 
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Bttring  No.  1 
L  ■  1.2S" 

D  »  1,25" 

C  ■  0,0005" 

m  a  8  X  10"^  lb, 


Bttrini  No,  2 
L  ■  1.25" 

0*2. 50" 

C  •  0.0010" 

m  B  12. 4  X  10“^  lb. -••c.Vtn- 


k  ealculAtad  from  Tig.  7 
■  1 . 43  X  10^  lb.  /in. 


Therefore 

r 

1 


2V  ^  -  675  -/MC. 


k  ciieulated  from  Fig,  7 

-  2.  12  X  10*  lb. /in, 

Therefore 

£  s  895  ^/eec. 
n 


Theee  compered  very  fevorebly  with  the  expertmental  value  of 
696  cyclei/eee. 

The  above  elarifiee  lome  of  the  obiervxtiona  and  chowa  comparatively 
good  agreement  between  theory  and  practice  In  the  caJculatlon  of  crltieal 
epeede  and  frequenelee  of  inatablllty. 
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X!i  conclusion:*  and  recommendations 

It  w&i  ihown  thAt  if  the  Re/noldi  EquAtion  it  integrAted  And  the  lllni 
fefeei  on  the  joumAl  computed  for  tero  tAngentiAl  velocity)  then  the  effects 
of  A  tAngentiAl  velocity  cAn  Also  be  obulned.  This  is  done  by  correcting 
the  rsdie!  velocity  (V^  <>  Cmc')  end  the  redlAl  And  tAngentiAl  force  compunenti 
(F^  And  F^)  using  the  fACtor  {1-ia').  ac  shown  in  EqoAtlons  (8  and  8a; 

Solutions  for  radlAl  and  tAngentiAl  forces  as  functions  of  e,  e*  And 
L/D  were  obtained  for  both  compressible  and  Incompressible  journal  bearings. 
Likewise,  the  derivatives  of  these  forces  with  respect  to  £  and  e*  were 
obtained. 

Calculations  were  also  made  of  the  forces  generated  in  externally 
pressurised  compressible  fluid  bearings.  It  was  shown  that  stiffness  in  this 
type  of  bearing  is  constant  for  at  least  small  eccentricity  ratios.  The  critical 
speed  is  independent  of  external  load.  It  was  further  shown  that  the  optimum 
bearing  from,  the  standpoint  of  maximum  stiffness  exists  at  L/D  ratios  around 

unity.  Miximum  stiffness  is  a  function  of  p  /p  ,  L/D,  and  Z  , 

sat 

Equations  of  motion  were  set  up  for  sn^AU  oscillations  about  a  position 
of  equilibrium  corresponding  to  some  particular  external,  steady  load.  These 
equA'nSr.3  contain  functions  f^,  f^  and  their  derivatives  with  respect  to  p  and 
e'.  These  differential  equations  are  linear  in  the  variable  V  ^'r  *•  y  and 
solutions  exist  in  which  time  enters  as  an  exponential.  Solutions  have  be-n 
obtained  for  both  compressible  i^nd  incompressible  cases  to  show  the  effect 
of  e  and  E  on  the  threshold  of  instability. 

Results  show  that  the  threshold  of  instability  Increases  with  Increase 
in  axtsrnal  load.  Rsftr  to  Flfures  ZZ  and  23. 
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TK«  ejLlen*ttle&  of  eHdefcl  tjMitd  dtiirvii  addltleaid  &tt4Btlon  tor 
e&iai  vhtn  tht  Ia  ail  diraotiQaa  aroto^  u  ■QuUtbrlum  pcint  li 

ttaa<iUAl. 

Tha  aqtutiosj  of  Urga  dlapUeamaBta  huTi  *lao  baaa  tat  in  tama 
of  e  and  a.  Mathod  for  obtaining  tha  path  of  Jeonal  eantar  for  tha  caaa 
of  a  maaalaaa  rotor  hava  batn  mttUnad. 

Tha  thraahold  for  Half  fraqaaney  Whirl  and  ethar  forma  of  dynamic 
difturbanea  hava  baan  traatad  in  aaetlona  vi  through  V’li  and  equations  are 
given  for  the  frequencies  of  disturbance. 

Wharayar  data  waa  ayailabla,  thaery  and  axparlaant  haa  baan  eomparad. 
Tha  raauiti  appaar  to  hart  good  agraamast.  Sayaral  of  tha  contrma ratal 
obiarrationi  hara  baan  elarlfiad. 

Tha  aabjaet  ii  by  no  maana  eemplataly  eorarad  and  It  ftUl  damanda 
eenaldarabla  atudy.  foraral  of  thaaa  atadiaa  art  anomaratad  balowi 
1)  Study  tha  baharler  of  rotora  for  imall  oaeillationa  at  ytry  low  accentriclty 

ratioa  and  whara  tha  angular  ineramaat  da  may  not  be  nagiscted. 

Z)  rind  raathodi  for  aolotion  of  tha  aquatloaa  of  nwtion  for  large  dla- 

plaeamanta.  Prafarably  try  to  obuin  eloaed  aelutloaa. 

3}  For  tha  caaa  of  large  diaplaeamanta  atudy  the  apaelal  caaaa  where 

a)  tha  accantrielty  ia  amall 

b)  tha  rotor  maaa  la  email  (mathod  of  aoletion  la  ootiinad  In  lactlon  Xi 

c)  tha  bearing  clearance  in  amall  eomparad  to  tha  abaft  banding 


daflaetlona. 
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4)  Study  the  above  three  pointi  in  caeei  where  the  rotor  or  bearing!  are 
not  eymmetrical. 

5)  Study  reethode  for  predicting  critical  epeeda  of  rotori  where  unequal 
•  tlffneiiea  in  two  direction!  exiat. 

6)  Conduct  experiment!  to  verify  all  the  above  point!. 

7}  Perform  experiment!  on  externally  preaaurixed  bearing!  to  determine 

how  eloiely  theory  agree!  with  practice  with  regard  to  conatant  itiXfneaa< 
optimum  L/D  ratio  and  critical  apeed  aa  function  of  eccentricity  ratio 
and  p  . 

8)  Study  pocketed  externally  preaiuriaed  journal  bearing!  under  the  aame 
condition. 

9)  Study  other  forma  of  inatability  which  have  been  reported  by  leveral 
isveatigaiora. 

10)  Eatabliah  chart!  and  table!  to  aid  englneera  in  the  deeign  of  liable 
bearing!. 
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IN’COMPRESSIBI-E  CASE  -  FULL  CIRCUL.*J^  SEARIKQ 
(Equal  Radial  and  Tangential  Shaft  Center  Velocitiee) 


TABIS  I 
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INCOMPRESSIBLE  CASE  -  FULL  CIRCULAR  BEARINQ 
{No  Tfcfigentitl  Shiit  Center  Velocity) 
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